Visualizing Spacetime Curvature via Gradient Flows I: Introduction 
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Traditional approaches to the study of the dynamics of spacetime curvature in a very real sense 
hide the intricacies of the nonlinear regime. Whether it be huge formulae, or mountains of numerical 
data, standard methods of presentation make little use of our remarkable skill, as humans, at 
pattern recognition. Here we introduce a new approach to the visualization of spacetime curvature. 
We examine the flows associated with the gradient fields of invariants derived from the spacetime. 
These flows reveal a remarkably rich structure, and ofi'er fresh insights even for well known analytical 
solutions to Einstein's equations. This paper serves as an overview and as an introduction to this 
approach. 

PACS numbers: 04.20.Cv, 04.20.Ha, 02.70.-c 



I. INTRODUCTION 

To quote from a recent monograph f]\: "It appears 
that it is much easier to find a new solution of Einstein's 
equations than it is to understand it". The present work 
involves the novel use of both computer algebra (for back- 
ground calculations) and numerical routines (for visual- 
ization) with the objective being the development of a 
fresh view of "curvature" . In this endeavor we are not 
alone. In 0] a visualization of spacetime curvature based 
on the electric and magnetic components of the Weyl ten- 
sor is being developed. In j3j] a potentially complemen- 
tary approach is also under development. Here we exam- 
ine something akin to a dynamical systems approach: we 
look at the "flows" associated with the gradient fields of 
invariants. We show that these flows reveal a remarkably 
rich structure and fresh insights even for well known ge- 
ometries. The present approach does not single out any 
particular observer, nor theory, but is designed to reflect 
deep intrinsic properties of the background geometry it- 
self. 



II. GRADIENT VECTOR FLOWS 

A. Invariants 

Let (Af, g) be a semi - Riemannian manifold of dimen- 
sion d|j where g is the metric tensor. The simplest con- 
struction of invariants involves scalars formed from met- 
ric contractions and partial derivatives to order p (scalars 
polynomial in the Riemann tensor). The number (Af) 
of algebraically (not functionally) independent scalars of 
this type (invariants not satisfying polynomial degenera- 
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for d > 2,p > 2 and so, for example, in spacetime the 
number of non-differential invariants (no derivatives of 
the Riemann tensor) is A/'(4, 2) = 14. This does not 
mean that a single set of 14 invariants will cover all pos- 
sibilities. The situation is rather more involved [^j and 
the issue of a minimal complete set remains unresolved 
even at rf = 4. Physical situations involve many simplifi- 
cations, and a consequent reduction in M . For example, 
it is well known that in the Ricci - flat case there are at 
most 4 invariants at order p = 2. Moreover, some invari- 
ants can be distinguished when their physical meaning 
is cleai0 and some can be used to construct Newtonian 
"analogues" as explained below. Whereas the mathemat- 
ics of polynomial invariants remains a very active area of 
research (see, for example, and references therein), 
as recently pointed out by Page 0, non-polynomial in- 
variants can be constructed even when all polynomial 
invariants vanish. 



B. Vector Fields 

Given (M, g) and a set of invariants X„ (we need not 
specify d nor p at this point nor even the nature of the 
In) consider the gradient flowf0 Q 



dir, 



(2) 



: 4 but higher dimensions 



For example, in Einstein's theory in 4-dimensions, without cos- 
mological constant and for a perfect fluid, Ti, = 87r(p — 3p), where 
'R. is the Ricci scalar, p the energy density and p the isotropic 
pressure. 

We refer to the choice " -|- " as a positive gradient flow and the 
choice " — " as a negative gradient flow. We find occasion to use 
both. 
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Note that k can be timehke, spacehke or null. Also note 
that even if a flow resides entirely in a subspace, calcula- 
tions must be carried out in the full space since invariants 
of the full space are not, of course, the same as invariants 
of subspaces. For notational simplicity we now consider 
n understood. By construction, we have 

V[afc/3] = ki^Vf^k^] = (3) 

and so on. In terms of the norm 

k^k" = V, (4) 

it follows from ([3]) that the acceleration associated with 

kiE 

k^Vpka. = iv„V. (5) 

From ([2]) it follows that the associated expansion is 

VaA:" = ± Va ( V"Z) = ±n T. (6) 

In terms of the projection tensor /i^ = (5J| ± k^-kp, 
borrowing the usual definition of the vorticity tensor 
ujaj3 = h1h^p\/[^ks], it follows immediately that 

= 0. (7) 

C. Special Symmetries 

For any scalar field I, the Lie derivative associated 
with a vector field is given by [s| 

/:eX = Cc.v"x = ea^"- (8) 

In certain circumstances, the scalar ^afc" can be directly 
related to the invariant I itself. For example, if the mani- 
fold admits a homothetic motion, C^gap = 2^ gap, where 
3" (= Va^"/d) is constant, it is known that for polyno- 
mial invariants [l^l 

C^I = K^I (9) 

where k is an integer characteristic of I (involving, for 
example, p and the number of discrete contractions used 
to make I). In the case that £, is Killing ($ = 0) we have 

Cafe" = 0. (10) 

That is, polynomial gradient flows are orthogonal to 
Killing flowsO This result can be used to clarify and 
extend the classical notions of TZ and T regions of space- 
time. This is explored in Appendix A. 



D. Special Invariants and Newtonian Analogues 

As mentioned above, in physical situations there is 
a reduction in Af due to simplifying assumptions. For 
example, in class B warped product spacetimes (which 
includes, for example, all spherical, plane and hyper- 
bolic spacetimes) it is known that A/'(4, 2) < 4 and that 
7V(4, 2) = 1 in the Ricci-flat case However, in addi- 
tion to simplifications leading to a decrease in the number 
of invariants that need to be considered, it has to be re- 
alized that not all invariants are equally useful. In this 
section, special invariants, that allow the construction of 
Newtonian analogues, are examined^ 

Decompose the Weyl tensor into its electric and mag- 
netic parts in the usual way: For a unit timelike vector 
w", Ea-y = Cap-ysu^u^ and H^j = Cap^su^u^ where C is 
dual to the Weyl tensor C. Since Hap has no Newtonian 
counterpart, consider here the simplest case, purely elec- 
tric spacetimes H^p = (e.g. all static spacetimes, shear- 
free and hypersurface orthogonal perfect fluid spacetimes 
etc.). Following for example Ellis f22l|, construct the New- 
tonian tidal tensor 

Eab = $,a,h - ^Vab^ $ (11) 

and define the associated invariant 

Sa6£;"'' = $,a,6*'"''-^ (□$)'. (12) 

We now construct two gradient fields 

k^ = ±V'^{EapE^P), l- = ±V\EabE-''). (13) 

We say that I is a Newtonian analogue (in no way any 
limit) of k if their associated phase portraits are "anal- 
ogous". By "analogous" we refer to, for example, the 
topology of the fiowQ To further clarify this construc- 
tion, let us recall that the tidal tensor in standard gen- 
eral relativity is given by [23'] Ta^ = TZaPjSU^u^ where 
again u" is a unit timelike vector and TZ is the Riemann 
tensor. In the Ricci - flat case then Ta-y — Ea-y and since 
CaPysC"'^'^^ cx EapE"!^ ^ the analogy constructed is 
strict since Cap-ysC"^''^ G I. Now, whereas V is con- 
structed in {M,g), is constructed in a Euclidean 3- 
space with metric rjab- Clearly the use of the term "anal- 
ogous" must account for this. The Curzon - Chazy solu- 
tion provides a detailed demonstration, and clarification, 
of these ideas [2^ . 

E. Singularities 

The gradient flows ([2]) can encounter a variety of excep- 
tional situations, the most common of which are critical 



It is essential that k is not normalized, since with jsj, normal- 
ization would restrict A; to a geodesic flow. 
^ In the special circumstance that k itself satisfies C^gaiB = 2^9^^ , 
it follows that k is concurrent Q, Vafc^ = ^ga/i, and so the 
associated streamlines are geodesic with V^V = 2^ka. 



^ A much simpler "analogue" is to simply construct the associated 
potential flow, the potential being the Newtonian potential. This 
is discussed below. 

Of course, the flow must have sufficient structure so as to be able 
to form an "analogy". 
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points were = 0. This is discussed below. Perhaps 
the simplest exceptional situation is the development of 
a caustic where the expansion of the flow diverges so that 
from ([5]) 



\ai\ 



— > oo. 



(14) 



In a recent influential work [26| , it has been claimed that 
attempts to explain cosmic acceleration, by way of obser- 
vations in locally inhomogeneous spacetimes, necessarily 
involves a local "weak singularity" at the origin. This 
argument is based on the observation that □ TZ neces- 
sarily diverges at the origin in these models. This point 
of view has been criticized (27| . In terms of a gradient 
flow, as introduced here, the correct way to view this 
controversy, is to observe that for I ~ TZ this "weak sin- 
gularity" is actually a caustic in the gradient flow of TZ, 
as follows from Since the models under consideration 
are simply dust, k represents the gradient flow of the en- 
ergy density and the "weak singularity" derives from the 
lack of sufficiently high differentiability of the energy den- 
sity at the origin. More substantial, and certainly more 
well known singularities, occur where an I diverges (so 
called "scalar polynomial singularities" in the simplest 
cases). In this circumstance, the "singularity" must be 
expanded to reveal the structure of the associated gra- 
dient flow about the "singularity" . In archetypical cases 
this reveals a very rich structure [1^ [1^ . 



III. DYNAMICAL SYSTEMS 

A. Contact With Physics 

Let us assume regular coordinates in a region of 
(M,g) and write the negative gradient flow in the form 



flows on (M, g) , yet most information readily available 
involves gradient flows on Euclidean n -space. This 
situation is addressed in Appendix B. In Appendix C we 
characterize fixed points by indices. 



IV. ELEMENTARY EXAMPLES 

A. The Robertson- Walker Spacetime 

Perhaps the most widely recognized spacetime is that 
of Robertson and Walker, 

rlr'^ 

ds^ = aitf{ —+r^dnj)-dt^. (16) 

1 — fcr^ 

Due to the isotropy, it is clear, a priori, that the associ- 
ated gradient "flows" can only be 1-dimensional. There 
can be no isolated critical points. The positive gradient 
flows are given by 



dt * 



so that 



(17) 



(18) 



The flows are timelike except at critical isotropic 3- 
surfaces for which dl/dt = 0. There are exactly two 
independent invariants at order p = 2, Ricci invariants, 
and allowing rational exponents of polynomial invariants, 
these invariants can be taken to be simply the energy den- 
sity p and isotropic pressure p. Turn - around Universes, 
those for which da/ dt\tQ = 0, deflne (irrespective of any 
cosmological constant) critical 3-surfaces t — to. 



dl 



(15) 



where ■ = d/dX and A is any natural parametrization of 
the curve a:" (A) with tangent k". It is clear that we are 
dealing with an autonomous dynamical system. In order 
to make an immediate contact with physics, we simply 
note here that the character of critical points associated 
with the system ()15|) for the Kerr metric is dependent 
on the value of the rotational parameter A [33| • There is 
"hidden" information intrinsic to the Kerr metric not ev- 
ident without the study of the associated gradient flows. 



B. Critical Points 

The phase portrait associated with is dominated 
by critical points, the deep study of which leads to Morse 
Theory |3l| . Gradient systems, in the large, have many 
diverse uses, and are the subject of much modern math- 
ematical study. Our considerations here involve gradient 



B. The Schwarzschild Spacetime 

Consider the Schwarzschild vacuum in dimensionless 
null coordinates {u,v,9,(j)) , 



ds^ = Am^ 



AC dudv 
{l + C) uv 



{l + Cfd^l 



(19) 



where m > is the effective gravitational mass, £ = 
C{uv) and C is the Lambert W function [sj] [35|. It 
remains convenient to define r = 2to(1 + C). It follows 
that ^" = {u, —V, 0, 0) is a Killing vector (and tangent to 
trajectories of constant r). We find 



{Am)' 



(2m-r),rV„e 



I r 

(20) 

and so along u — Qorv — Q{r — 2m), ^ is tangent 
to a radial null geodesic. Further, we find that Vq,^" = 
i[cx'^ p^-y] = so that regions for which r > 2m are static; 
all of which is, of course, well known. 



4 



For (I19p . there is only one independent invariant of 
order p — 2, and we can take this to be the square of the 
Weyl tensor 



1 = 



4to4(i + c)6 



(21) 



It fohows that the tangents to the negative gradient flow 
are given by cx —m{u, v, 0, 0)/r^ (tangent to trajecto- 
ries of constant "Schwarzschild t" ) , 

V oc — (r - 2m), k^V^k^ cx jg-^^ (u, w, 0, 0), 

(22) 

and Vak'^ cx m^(5r — 12m)/r^. Again, along u — 
OT V = 0, k is tangent to a radial null geodesic. The 
TZ. region corresponds to r > 2m whereas the T region 
corresponds to r < 2m, exactly as expected. Both ^ and 
k vanish identically at the isolated critical point V, the 
bifurcation 2-sphere u = v = [3^ . 

Let us now be more explicit about the flow with the 
aide of Appendices B and C. We find 



where 



( 3 



V4mV {l + C)"' 



In particular, we find 



and 



dG 



3m" 



4m^ 



> 



> 



(23) 



(24) 



(25) 



(26) 



and so the isolated critical point is a stable node and the 
associated index is -|-l|f| In the neighborhood of uv = 0, 
it follows from (l24l) that -0 ^ — (3/4m^)^ and so from 
(|23p it follows that we have a simple linear product flow 
in the neighborhood of uv = 0. Such fiows have a natural 
interpretation in terms of gradient flows |33l ■ 

The phase portrait in this case is very simple and the 
associated Newtonian potential flow is easily constructed. 
Setting the center of symmetry at the bifurcation two- 
sphere, and using Cartesian coordinates, it follows that 
the associated Newtonian potential flow is that of a fluid 
sphere of uniform density. In situations with more sub- 
stantial phase portraits, one can go beyond associated 
potential fiows and construct Newtonian analogues as ex- 
plained above. 



As mentioned in Appendix A, the introduction of a cosmologi- 
cal constant A alters the topology of the conformal boundary of 
spacetime. In particular, the Schwarzschild de - Sitter case has 
two bifurcate 2 - spheres and a consequent index of +2. 



C. Weyl Invariants and Spherically Symmetric 
Spacetime 

As mentioned above, in class B warped product space- 
times it is known that A/'(4, 2) < 4 and that A/'(4, 2) = 1 
in the Ricci-fiat case [2l|. Here we examine Weyl invari- 
ants within spherical symmetry. (In the notation of f2l| , 
we need only consider the first Weyl invariant wi (equiv- 
alent to the square of the Weyl tensor, divided by 8).) 
Now locally^ every spherically symmetric spacetime can 
be written in the form 



-2/(u, v)dudv + r{u, v)^dn2 



(27) 



It follows that the negative gradient fiow of wi is given 

by 



with 



where 



and 



A:" = -A(u, v){fi{u, v), d{u, v),Q, 0) 



V cx -A^ffiS 



wi cx AV^ 



3/8r5 ■ 



(28) 



(29) 



(30) 



(31) 



^ , fi and S involve derivatives of / and r up to order p — 
2. The explicit forms are not needed here. Rather, we 
note that there exist two distinct types of critical points 
associated with the gradient flows of wi. From (|28l) and 
(1311) it follows that the critical points can be characterized 
by the local conditions 



^- = 0, ^ = 5 = 0. 



(32) 



In the Robertson- Walker case, of course, 5* = is a global 
condition. The Schwarzschild spacetime demonstrates a 
case where ^ locally. Rather, the bifurcation 2- 
sphere corresponds to an isolated critical point charac- 
terized by fi\-p = S\-p = 0. We call such critical points 
locally anisotropic. 

In order to explicitly demonstrate examples where = 
locally, consider, for example, static perfect fiuids in 
comoving coordinates. 



ds' 



dr^ 



1 



2m(r) 



r^dnj 



(33) 



Given $ sufficiently smooth, and subject to boundary 
conditions, one can generate the effective gravitational 
mass m that gives rise to a perfect fluid [s^l . For example, 
setting 



1 



7Vln(l + — ) 



(34) 
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where N is an integer > 1 and a is a constant > 0, an 
infinite number of exact perfect fluid solutions follow. For 
all such models we have 

= V = □ wi = (35) 

at the center of symmetry r = 0. Moreover, k" = 
(0, 0, 0, 0) there and so we have a critical 3-surface at 
the origin based on the Weyl invariant. That is, the so- 
lutions are locally isotropic about the origin for all t. We 
call such a collection of critical points locally isotropic. 

V. DISCUSSION 

This paper serves as an introduction to the study of 
gradient flows of invariants as a means to visualize cur- 
vature. No particular observer, or theory, is fundamental 
to this approach, but the intention is to apply the tech- 
niques to spacetimes associated with Einstein's theory of 
gravity. It has been shown that gradient flows of invari- 
ants polynomial in the Riemann tensor are necessarily or- 
thogonal to Killing flows, should they exist. This allows 
for a rigorous deflnition and generalization of the classi- 
cal notions of TZ and T regions of spacetime. For purely 
electric spacetimes one can construct strict Newtonian 
analogues, based on the topology of flows; for the first 
Weyl (tidal) invariant in spacetime, and for the tidal in- 
variant in Newtonian theory. This is discussed at length 
elsewhere ^25]. 

A central feature of many spacetimes is the develop- 
ment of "singularities" , a concept which remains open 



to much refinement. It has long been known [38l| that 
naked spherically symmetric singularities are "massless" , 
a fact clearly reflected in the associated gradient flows 
[28*1. This invites a reexamination of the very notion of 
"singularity" . 

From the point of view of dynamical systems, gradi- 
ent flows are simple in the sense that critical points are 
the only possible limit sets of the flow. We have trans- 
lated standard results in into the setting of a semi- 
Riemannian manifold and reviewed the importance of the 
index associated with each critical point, reflecting on the 
role of the cosmological constant as regards its role on the 
fundamental index of a spacetime. 

We have reviewed examples, restricted here to spheri- 
cal symmetry. A rather complete analysis of the Kruskal 
- Szekeres vacuum was given, interpreting the associated 
bifurcate 2-sphere as the isolated critical point of the so- 
lution. More generally, isolated critical points of the Weyl 
invariant within spherical symmetry were distinguished 
as locally isotropic or anisotropic, with explicit examples 
given for each. 

Finally, we note that even at the level of complexity of- 
fered by the Kerr spacetime, a remarkably rich structure 
is revealed by the associated gradient field [30]. 
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Appendix A: TZ and T regions 

The gradient of the areal radius of an arbitrary spheri- 
cally symmetric field can be spacelike, timelike or null. A 
careful distinction of these possibilities is important and, 
for example, forms an essential element of an y co mplete 
proof of the Birkhoff theorem (for example [il|). The 
possibility of spacelike and timelike gradients was stud- 
ied extensively in the Russian literature (and labeled "7?." 
and "7^' regions respectively) some forty years ago jl^j 
and yet there appears to be no readily available exten- 
sion of these ideas. It is the purpose of this Appendix to 
explore this using gradient flows. 

The original distinction of TZ and T regions of space- 
time was restricted to spherical symmetry and given by 
Novikov fl3\ (see also the work of Ruban fl^ ) by way of 
a coordinate construction J^3] . This construction can be 
made invariant as follows: Consider a spherically sym- 
metric spacetime [l3| 

ds^ = dsl + R^dnj (Al) 

where c?ri| is the metric of a unit 2-sphere {d9^ -\- 
sin^(0)d(/)^) and R = R{x^,x^) (the areal radius) where 
the coordinates on the Lorentzian 2-space S are labeled 
as and x^. (No specific choice of coordinates on E need 
be made.) Writing the effective gravitational mass (some- 
times called the Misner-Sharp energy) as m = m{x^ , x'^) 
it follows from Novikov's coordinate definitions that 

R>2m, R=2m, R< 2m (A2) 

in an TZ region, on the boundary and in a T region re- 
spectively [31 ■ In more modern notation [l7j . an TZ re- 
gion is untrapped, a T region trapped and the boundary 
marginal. Whereas these distinctions are fundamental, 
they are restricted to strict spherical symmetry. The nat- 
ural extension of these ideas beyond spherical symmetry 
involves finding the apparent horizon(s) [ll| . a key ele- 
ment of numerical relativity [18]. 

In order to bring gradient flows into the picture, let us 
note that a stationary region of spacetime admits a time- 
like Killing congruence. Since every non-zero 4-vector 
orthogonal to a timelike 4-vector must be spacelike, it 
follows from (jTU]) that any gradient flow is necessarily 
spacelike in a stationary region. From (U) we have V > 
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throughout an TZ region, a fact that we can use to define 
such a region. Now let us assume that there is a region 
in which the gradient flow is timehke so that V < 0. 
(A Kilhng congruence in this region, should one exist, is 
then necessarily spacelike.) We can use V = to define 
a T region. Boundary regions are then naturally defined 
by V = 0. The foregoing classification is independent 
of symmetries and we adopt this classification in general: 

V > throughout an TZ region where k is spacelike, V < 
throughout an 7" region where k is timelike, and V = 
on a boundary where k is null or the zero vector. It is 
already known [l^ , for example, that for the Kerr metric 

V > for the rotational parameter A > \ with equality 
holding only in the degenerate case A ~ 1. That is, de- 
generate and naked Kerr metrics do not admit T regions. 
Such statements are, however, strongly dependent on the 
sign of the cosmological constant (A). If A > then con- 
formal infinity is spacelike 20] and so one would expect 
that the degenerate and naked cases would then admit 
no TZ regions. 



Appendix B: The Euclidean Counterpart 

In this Appendix we wish to convert the covariant 
statement 



dl 



into the non-covariant but common form 

dF 



ax" 



(Bl) 



for the purpose of exploiting many well-known properties 
of the system (|B1[) . Here we restrict our considerations to 
2-dimensional flows, but we note that higher dimensional 
cases follow analogously. From ((T5|) we have 



dl 



We now simply write out the associated system 

r)F f)F 



(B2) 



(B3) 



where the identifications for F are obvious, and the last 
equivalence is for notational convenience. Note that we 
never need solve for F. Rather, we proceed in a purely 
mechanical way to construct the Hessian and its discrim- 
inant D. We need not write these out explicitly here. 



if the discriminant vanishes, V\-p = 0. Otherwise P is 
non-degenerate, that is, a Morse critical point. We are 
now in a position to summarize the situation 32| : Morse 
critical points of a gradient flow are hyperbolic. Further, 
for isolated Morse critical points we have: 

• "P is an unstable node if and only if 2?|p > and 

• P is an asymptotically stable node if and only if 
V\-p > and dG/dx^)\v > 0. 

• P is a saddle point if and only if Pjp < 0. 

It is to be noted that gradient flows are simple flows in 
the sense that critical points are the only possible limit 
sets. However, since we must also deal with singularities 
in the invariants I, we must also deal with sources and 
sinks where one or more of the diverge. 

Now whereas we do not have to solve for F, we can 
attempt to do so. In particular, we can write 



(B4) 



A critical point V, where x^ 



0, is degenerate 



If we take $ to be the Newtonian potential, we refer to 
the flow (|B4|) as the associated Newtonian potential flow. 



Appendix C: Index Theory 

Again here we restrict our considerations to 2- 
dimensional flows, higher dimensional cases follow anal- 
ogously. We wish to characterize the topology of the 
2-surface associated with the flow [s^. We note the fol- 
lowing regarding the Poincare index: 

• The index of a node is -1-1. 

• The index of a source or sinl0 is -1-1. 

• The index of a hyperbolic saddle point is —1. 

• The index of a closed curve containing fixed points 
is equal to the sum of the indices of the fixed points 
within. 

The importance of calculating theses indices follows from 
the Poincare - Hopf theorem which we state here in the 
following simple form: Suppose the fiow is on a 2-surface 
E. Calculate the sum of the indices of all fixed points 
on S with a suitable curve 7 which contains all isolated 
fixed points. Then this sum is the Euler Characteristic 
of E within 7. 



These are included so as to be able to cover singular fixed points 



where one or more of the x°' diverge. 



